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ABSTRACT 

SDSSJ125733.63+542850.5 (hereafter SDSSJ1257+5428) is a compact white dwarf binary from 
the Sloan Digital Sky Survey that exhibits high-amplitude radial velocity variations on a period 
of 4.56 hours. While an initial analysis suggested the presence of a neutron star or black-hole binary 
companion, a follow-up study concluded that the spectrum was better understood as a combination 
of two white dwarfs. Here we present optical spectroscopy and ultraviolet fluxes which directly reveal 
the presence of the second white dwarf in the system. SDSSJ1257+5428's spectrum is a compos- 
ite, dominated by the narrow-lined spectrum from a cool, low gravity white dwarf (T c g ~ 6300 K, 
logg = 5 to 6.6) with broad wings from a hotter, high-mass white dwarf companion (11,000 to 
14, 000 K; ~ 1M©). The high-mass white dwarf has unusual line profiles which lack the narrow cen- 
tral core to Ha that is usually seen in white dwarfs. This is consistent with rapid rotation with 
vsini = 500 to 1750 km s -1 , although other broadening mechanisms such as magnetic fields, pulsa- 
tions or a helium-rich atmosphere could also be contributory factors. The cool component is a puzzle 
since no evolutionary model matches its combination of low gravity and temperature. Within the con- 
straints set by our data, SDSSJ1257+5428 could have a total mass greater than the Chandrasekhar 
limit and thus be a potential Type la supernova progenitor. However, SDSSJ1257+5428's unusually 
low mass ratio q sa 0.2 suggests that it is more likely that it will evolve into an accreting double white 
dwarf (AM CVn star). 

Subject headings: supernovae: general - white dwarfs - accretion - binaries: close 



1. INTRODUCTION 

Close pairs of white dwarfs with combined masses 
greater than the Chandrasekhar limit have long 
been discuss ed as potential proge nitors of T ype l a 
supernovae (|Iben fc Tutukovi 119841: iWebbinld [1 984). 
Howe ve r despite se a rches, ([Robinson & Shaftcr 
19871: IBragadia et al.l 11990: IMarsh et al.l ll995T 
Napiwotzki et alj I2003D no secure examples of dou- 
ble white dwarfs both massive enough and short-period 
enough to merge within a Hubble time have been found. 
Only about 1 in a 1000 white dwarfs ar e required to be 
in such systems to match Type la rates (jNelemans et alJ 
I2001bfl . ' Since only of order 1000 white dwarfs have 
been searched for binarity to date, the deficit is not yet 
significant, however it continues to be worth searching 
for more such systems. In an effort to do so we exam- 
ined the spectra of DA white dwarfs (those showing 



spectra with hydrogen absorpt ion only) from the SDSS 
survey (jEisenstein et al.l 120061 ). looking for objects of 
discrepant radial velocity. One star, SDSSJ1257+5428 
was an obvious outlier with a mean radial velocity 
of — 300 km s . In this paper we present follow-up 
spectroscopy to elucidate the nature of this object. 

S DSSJ1257+5428 was the subject of a similar study 
by iBadenes et"aH (|2009f ). They found that it was a 
binary, measured a radial velocity semi- amplitude of 
K\ = 323 ± 6 km s _1 on a period of 4.56 hours and fitted 
their spectra with a white dwarf of temperature ~ 9000 K 
and high mass, ~ 0.9 M©. The period, radial velocity 
amplitude and white dwarf mass give a minimum mass 
for the companion of 1.6 Mq, suggesting tha t it is either 
a neutron star or black-hole. IBadenes et al.1 (|2009f ) esti- 
mated a distance of 48 pc for SDSSJ1257+5428, implying 
t hat such systems may be rather common. 

IBadenes et al.l (|2009f ) recognised that their spectral fits 
were problematic, and in particular failed to fit the nar- 



2 




-500 500 -500 500 -500 500 -500 500 



-2000 2000 -2000 2000 -2000 2000 -2000 2000 



Velocity (km s ') 



Velocity (km s ') 



Fig. 1. — The panels show phase-folded trailed spectra of the first four Balmer lines of SDSSJ1257+5428, Ha— 5, left to right. The 
left-hand panel shows the raw data; the right-hand panel, which shows the data after removal of the primary's motion and mean spectrum, 
reveals anti-phased features from the massive secondary white dwarf. (Note the change of horizontal scale between the two panels.) 

row cores of the Balmer lines. Thus while they favored 
a neutron star or a black-hole for the companion, they 
could not entirely eliminate the possibility that it w as an- 
other white dwarf. iKulkarni fc van Kerkwiikl ()2010[ ) took 
three high signal-to-noise spectra which showed asym- 
metries suggesting exactly this; they showed that their 
spectra could be fit by a combination of a cool, low 
mass white dwarf (6250 ± 250 K, 0.15 ± 0.05 M ) plus a 
hotter, massive white dwarf companion (13,000 ± 800 K, 
0.92 ± 0.13 M, 



IKulkarni fc van Kerkwiikl (|2010D 's paper appeared 
shortly after the original submission of our work. Al- 
though both our papers agreed upon the basic double 
white dwarf nature of SDSSJ1257+5428, there were sig- 
nificant differences of detail, with inconsistencies in both 
masses and temperatures. In an effort to understand 
these, we have since carried out additional fits to our 
data, which we present here along with our original ap- 
proach. In addition, an improved Swift calibration has 
given us more confidence in the modelling of the UV- 
optical spectral energ y distribution. Our new results 
agree more closely with Kulkarn i fc van Kerkwiikl (j2010f ) 
than our original analysis, but also reveal uncertainties in 
the system properties which make it impossible to estab- 
lish securely such fundamental properties as whether the 
system is super-Chandrasekhar or not. We begin with a 
description of our observational material. 

2. OBSERVATIONS AND REDUCTION 

We obtained 140 spectra in the B and R bands covering 
the Balmer series on the William Herschel Telescope in 
La Palma using the ISIS double-beam spectrograph. The 
first 22 spectra were taken in service mode on the night of 
April 20, 2008, while the final 118 were acquired during 
a four night run spanning April 29 to May 2, 2009. The 
data were debiassed, flat- fielded and optimally extracted 
(|Hornell 19861: lMarshll989l) . The wavelength scale for each 
spectrum was linearly interpolated in time from pairs of 
bracketing arc calibration spectra. In 2008, the night was 
clear with seeing 0.8 to 1.6 arcsec, however the Moon was 
full. In 2009, the Moon was only 40% illuminated, but 
conditions were poorer, and we lost a total of one night 
to clouds while the seeing varied from 1.0 to 2.5 arcsec. 

In addition to the optical spectroscopy, we obtained 
ultraviolet flux measurements of SDSSJ1257+5428 with 
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Fig. 2.— Radial velocities of Ha from 2008 and 2009 folded on 
the ephemeris of this paper. 

the Swift satellite and accessed additional archival data 
from 7 epochs of observations covering the interval Aug 
24 2009 - Jan 13 2010. Our summed images had 
effective exposure times ranging from 818 to 5608 s. 
SDSSJ1257+5428 was well detected in all four U V filters 
provided by Swift /UVOT (jRoming et all 120051 ). Aper- 
ture photometry with a 5" extraction region and an an- 
nular sky region was performed using the SWIFT data 
analysis package V3.7 within the HEASoft V6.10 release. 
Associated zero-point calibration files and filter effective 
area curves (Jan 2011 release) were used to calculate the 
fluxes for each filter. 

3. RESULTS 

In all that follows we will refer to the component most 
obvious in the spectra as the primary star, and its com- 
panion as the secondary star. This means that the pri- 
mary star is low mass and cool, while the secondary star 
is high mass and hot. 

3.1. The primary star 

Orbital motion of large amplitude was immediately ap- 
parent in the spectra. The left panel of Fig. [1] shows 
phase-folded trailed spectra of the first four Balmer lines 
displaying high-amplitude orbital motion and an appar- 
ently single-lined binary star. 
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Fig. 3.— The radial velocity semi-amplitude of SDSSJ1257+5428 
from Ha to HIO showing the varying influence of the companion 
(which always reduces the amplitude) across the Balmer series. 

We fitted the radial velocities (Fig. [5]) of the Balmer 
lines using combinations of 3 Gaussians for each line, 
with full width half maxima (FWHM) and depths op- 
timised to minimise \ 2 - Given the large radial veloc- 
ity amplitude of SDSSJ1257+5428, the four night run 
in 2009 was sufficient to extrapolate back to the service 
night of a year earlier to find a unique alias, leading to 
the following ephemeris (on a UTC timescale corrected 
for light- travel time to the heliocentre) : 

HJD = 2454846.17470(8) + 0.18979154(9)£, (1) 

where the time of zero phase corresponds to the time 
when the primary star is closest to Earth. Wh i le our 
period is consistent with that of IBadenes et al.l (2009) 
(P = 0.18979(3) days), our zero phase is not, since on 
our ephemeris their zero phase occurs at E = 115.619 ± 
0.004; this is the result of an error of 0.5 days in the 
calculation of Julian Days together with no allowance 
for light travel time in IBadenes et al.l ([2009D (C.Badenes, 
priv. comm.); we have in addition confirmed that our 
ephemeris correctly predicts the velocities of the SDSS 
spectra which were taken in 2003. 

Having established a unique period alias, from now 
on we use only the 2009 data because of differences in 
the instrumental setups between 2008 and 2009 which 
make consistent continuum fitting difficult. Sinusoidal 
fits were made to the radial velocities measured for each 
Balmer line from Ha to H10. The semi-amplitudes were 
seen to vary from 324.1 ± 0.8 km s -1 for Ha to 300.6 ± 
2.6 kms -1 for He, rising again in the higher Balmer lines 
(Fig. [3]). The variation from line to line is a symptom 
of the presence of a second broad-lined white dwarf. As 
indicated earlier, and as will be seen later, the secondary 
white dwarf is hotter and more massive than the primary. 
Therefore, owing to its temperature on the one hand, and 
its high gravity on the other, the secondary is expected 
to contribute relatively weakly at Ha and in the higher 
Balmer series, but more strongly in intermediate series 
lines. This is exactly what is seen in Fig. [31 with the 
radial velocity semi-amplitude acting as a proxy for the 
relative line strength from the two components. 

3.2. The secondary star 

To search for more direct signs of a companion we 
first fitted and normalised the continua of the spectra 
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Fig. 4. — The average _B-band spectrum (top) and results of dis- 
entangling the spectra of SDSSJ1257+5428 (gray lines) along with 
model atmospheres (primary: T e ff = 6900 K, logg = 6.75; sec- 
ondary: T cff = 10300 K, logg = 8.70) for a light ratio of 0.7 (sec- 
tion !3.4l l and usini = 1500 km s — 1 (black lines). For disentangling, 
the radial-velocity semi-amplitudes were fixed at K\ = 330kms — 1 
and K2 = 100 km s" 1 , but the results were insensitive to the pre- 
cise value of K2- 
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Fig. 5. — The disentangled Ha spectra with the secondary star 
component plotted at the top. As with Fig. [4] the continuum levels 
are indeterminate. 

between the Balmer lines. Next we shifted out the pri- 
mary star's orbital motion and subtracted the average 
of the resulting spectra. If only one component con- 
tributes significantly, this procedure should leave only 
noise. Any contribution from a secondary star will cause 
correlations amongst the residuals. The right-hand panel 
of Fig. [1] shows such residuals which moreover are anti- 
phased compared to the left-hand panel. This is the 
signature of a second contributing white dwarf, albeit 
distorted by the subtraction process. 

The spectra of two stars within a binary can 
also be recovered using a method applied successfully 
to main-sequence bin aries known as "disentangling" 
(|Simon fc Sturm|[l994D . Although the broad absorption 
lines of white dwarfs make it harder to establish the reli- 
able normalisation required for disentangling, we tried it 
out to see if any separation of the two spectra was possi- 
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ble. Its appeal is that it requires no a priori assumptions 
concerning the spectra of the two components, although 
it cannot return the individual continua, and takes as 
input continuum-normalised spectra. The results for the 
blue spectra, plotted in Figs 0] and over-plotted with 
model atmospheres, show a pair of low- and high-gravity 
white dwarf spectra. The high-gravity secondary con- 
tributes the b r oad l ine wings and it is thus clear why 
iBadenes et alj (|2009T ) over-estimated the mass of the pri- 
mary and had difficulties obtaining good spectral fits. 
The process of disentangling leaves the continuum levels 
undetermined; we discuss the model atmosphere fits fur- 
ther in section 15^41 The disentangled Ha profiles (Fig. [5} 
are remarkable for an absence of any sharp line core in 
the high gravity spectrum. Ha lines in white dwarfs 
of this temperature usually display a sharp dip at the 
line ce ntre 2 or 3A in width as the result of non-LTE 
effects dGreenstein fc Peterson! Il973; Ko e ster fe Herrerol 
119881: IHeber et alj|1997l ), but Fig. [51 shows no sign at all 
of this. The figure also makes it clear why Ha relatively 
faithfully traces the motion of the primary. 

3.3. Mass constraints 

We were unable to obtain a reliable direct measure- 
ment of the semi-amplitude of the secondary star K% 
from our data, a consequence, we believe, of its high 
gravity and the absence of any sharp line cores, combined 
with the difficulty of defining reliable continua. However, 
when trying to obtain simultaneous fits to both stars, we 
found that the extra freedom of a secondary component 
always allowed K\ to increase above the directly-fitted 
values to a value consistently close to 330 ± 2kms _1 
which we take to be the true value of K\. The di- 
rectly measured values are biassed towards lower val- 
ues by the contribution from the secondary, although 
as Fig. [3] demonstrates, the effect is relatively slight in 
Ha and in the higher B aimer series where the secondary 
star's high gravity ensures that it contributes little. The 
value of Ki, along with the presence of a second white 
dwarf leads to the mass constraints plotted in Fig. [8] 

3.4. Spectroscopic fits 

Our WHT data reveal that the narrow Balmer lines 
in SDSSJ1257+5428 extend up to H12 (Fig. g]), im- 
plying unambiguously that the cooler white dwarf (pri- 
mary) is of an unusually low surface gravity. In 
the original submission of this paper, we modelled 
our WHT spectr oscop y of SDSSJ1257+5428 as we de- 
scribe in Sect. l3~4Tl finding that SDSSJ1257+5428 
is composed of a very low-mass, cool white dwarf 
plus a hotter and rather massive white dwarf that is 
rapidly rotating. While our pa per was under review, 
iKulkarni fc van Kerkwiikl ((2010h published an indepen- 
dent study of SDSSJ1257+5428, coming to broadly sim- 
ilar conclusions, but with noticeable differences in the 
detailed white dwarf parameters. In the light of this 
development, we re-visited the modelling of the WHT 
spectroscopy and the ultraviolet (UV) /optical spectral 
energy distribution (SED), adopting three different ap- 
proaches which we now detail. 
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Fig. 6. — Constraints on the effective temperatures and surface 
gravities of the two white dw arfs in SDSSJ1257+5428. Fits to 
the disentangled spectra (Sect. 1X4, 111 were carried out with light 
ratios 0.7—1.3 in steps of 0.1, and secondary star rotation rates 
vs'mi = 500 — 1750 kms" 1 in steps of 250 kms -1 , and are shown 
as open circles. The best-fit paramete rs from fitting the phase- 
resolved 2009 WHT spectroscopy (Sect. l3T4.2H are shown as filled 
circles, their averages and standard deviations are indicated by 
the red error bars (T ljeff = 6341 ± 139 K, log pi = 5.40 ± 0.37, 
r 2iCff = 11111 ± 704 K, and logg 2 = 9.13 ± 0.14). Fits to 
the S DSS/U VOT broad-band photometry are shown as triangles 
(Sect. | 3.4.3C. The blue error ba rs show the best-fit parameters of 
Kulkarni & van Kcrkwiik (2010), and the dashed line indicates the 
occurrence of maximum Balmer-line strength. 



TABLE 1 

Summary of effective temperature and gravity 
measurements for sdssj1257+5428. 



Source 


T e g : i 


log gi 


Teff, 2 


log 32 




[K] 


[cms -2 ] 


[K] 


[cms -2 ] 


Disentangling 


7200±350 


6.9±0.1 


9800±1000 


9.0±0.4 


Phase-resolved 


6340±140 


5.4±0.4 


11100±710 


9.1±0.2 


UV-optical SED 


6620±200 


6.1±0.1 a 


12400±710 


8.6±0.2 


KvK2010 


6250±250 


6.0±0.3 a 


13000±800 


8.5±0.2 



3.4.1. 



spectra 



1 Fits limited to logg > 6.0. 

We modelled the two individual spectral components 
resulting from t he disentangling of Fig. H] using the fit- 
ting routine of iRebassa-Mansergas et alj (|2007f ) . The 
gri d of DA model spe ctra which was calculated as 
in iKoester et alj (|2005() . include s the updated Stark 
broad ening of the Balmer lines of iTremblay fc Bergeron! 
(2009), and spans 5 < logg < 9.5 in steps of 0.25 and 131 
temperatures sampling the range 6000 < T e g < 100000 K 
in an optimised way. 

Given that the relative flux ratio of the disentangled 
components is indeterminate, we carried out fits adopt- 
ing a range of "light ratios", I = 0.7 — 1.3 (in steps of 
0.1), defined such that the secondary star contributes I 
while the primary contributes 2 — I to the observed light. 
Changing the light ratio modifies the ratio of the Balmer 
line equivalent widths between the two components, and 
light ratios outside the adopted range resulted in unphys- 
ical line strengths for one or both components. In addi- 
tion, while spectral disentangling adopts a single value 
for the light ratio across the entire spectral range that is 
analysed, the true light ratio of the two white dwarfs in 
SDSSJ1257+5428 is a function of the wavelengt h. Thi s 
may introduce systematic uncertainties (see Sect. l3.4.l2|) . 
although one would hope that the range of light ratios 
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Fig. 7. — Composite model fits (red) to the 2009 WHT spectroscopy of SDSSJ1257+5428, averaged in 30 orbital phase bins (black). 
Each individual spectrum was fitted independently, resulting in 30 solutions for T\ jE ff, loggi, T2 jC ff, log 32 (see Fig. [6}, and the relative flux 
contribution of both white dwarfs. 



adopted would reveal such systematica through variation 
in the best fit parameters. Conditions during our ob- 
servations precluded any reliable flux calibration of our 
spectra and we use no information from the continuum. 
For each assumed light ratio I, the model fit to the pri- 
mary results in a best-fit effective temperature and sur- 
face gravity, an example for I = 0.7 is shown in Fig. 01 
Figure[S] illustrates that regardless of I , the parameters 
of the primary are constrained within a narrow range, 
T cfU = 7200 ± 350 K and log 51 = 6.85 ± 0.10. This 
temperature and gravity impli es a very low mass for the 
primary, M x = 0.20±0.05 M (jPanei et al.ll2007D . These 
fit parameter s and others discussed in this pap er, includ- 
ing those of iKulkarni fc van Kerkwijkl (|2010l ) listed as 
KvK2010, are summarised in Table [Tl 

The broad Balmer lines of the secondary (particularly 
Ha) suggest rapid rotation. Assuming this, we found 
that acceptable fits were only obtained for projected 
equatorial rotation speeds in the range v sin i = 500 - 
1750 km s -1 . The distributions of the best-fit parame- 
ters in the T e g — log 5 plane are again illustrated in Fig. [SI 
with higher values of v sin i requiring lower log g to repro- 
duce the observed broad line profiles. These span a range 
T eff , 2 = 9800 ± 1000 K and log 52 = 9.0 ±0.4, correspond- 
ing to M2 = 1.2 ± 0.2 Mq, consistent with the kinematic 
constraints of Fig. [5J These parameters broad l y agre e 
with the findings of IKulkarni fc van Kerkwiikl (12010) . 
however, they found an even lower gravity for the low- 
mass white dwarf primary (loggi = 6.0 ± 0.3), and 
a higher temperature for the more massive secondary 
(T effi2 = 13000 ± 800 K). We will see below, as is 



also clear from Table [TJ th at the Swift UVOT fluxes 
of SDSSJ1257+5428 (section 13X3} indeed suggest that 
T e &,2 from modelling the disentangled spectra is too low. 

3.4.2. Phase-binned flux spectra 

An assumption inherent to the spectral disentangling 
is a constant light ratio, i.e. relative flux contribution 
of both components, across the wavelength range under 
analysis. Our application of the method over the en- 
tire wavelength range of the blue WHT spectra violates 
this assumption, which will result in systematic errors 
in the relative strengths of the individual Balmer lines, 
although as we remarked, our use of a range of light ra- 
tios might be expected to transfer any resulting variation 
into the final uncertainties. However, as an alternative 
we decided to carry out more direct fits to the WHT spec- 
tra. We binned the 118 individual spectra into 30 phase 
bins, and fitted each of them separately with the sum 
of two white dwarf model spectra drawn from the same 
model grids as in Sect. l3.4.T1 Observed and model spectra 
were normalised prior to the fits. The fits were repeated 
for a range of rotational velocities for the more massive 
white dwarf, 500 < wsini < 1750 km s -1 , in steps of 
250 km s _1 . 

Figure [7] shows the independent best-fits to the 30 
phase-binned spectra for wsini = 1000 km/s. The distri- 
butions of the 30 individual best-fit white dwarf param- 
eters are shown in Fig. [51 Each of the 30 independent 
fits leads to a pair of points in this figure, one in each 
of the two separate groups of points. Averaging these 
groups leads to T cfU = 6340±140K, log 51 = 5.40±0.37, 
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Fig. 8. — Constraints upon the masses of the two stars in 
SDSSJ1257+5428. The solid curved line is the lower limit upon M 2 
given K\ = 330kms — . The value of M2 is bounded from above by 
the Chandrasekhar limit (horizontal dashed line), Mch = 1-37 Mg, 
while a conservative minimum white dwarf mass of 0.1 Mg places 
a lower limit upon Mi . These three constraints lead to the shaded 
triangular region. The dotted-line box outlines the more arguable 
region of parameter space favored by the spectroscopic and evolu- 
tionary arguments outlined in section [4.1l 



T eff , 2 = 11100 ± 710 K, and log .92 = 9.13 ± 0.14. Higher 
(lower) values for v sin i shifts both white dwarfs to lower 
(higher) T e g and logg, by a few hundred degrees and 0.1- 
0.2 dex, with no significant difference in the quality of the 
fits. 

Compared to the fit of the disentangled spectra 
(Scct. l3"Al"T) . the two white dwarfs move further apart 
in the (T e ff,\ogq) plane , som ewhat in the sense of 
iKulkarni fc van Kerkwijkl ((2010). Interestingly, log 51 = 
5.40 ± 0.37 from fitting the ph ase-binned spectra is even 
lower than that reported by IKulkarni fc van Kerkwijkl 
(2010), log gi — 6.0 ± 0.3. However, their best-fit value 
for log gi is on the boundary of their model grid, and thus 
may have been limited by the range of their models. This 
exceedingly low gravity by white dwarf standards raises 
evolutionary problems as we elucidate in section |4. II 

The best-fit parameters show considerable scatter in 
the gravity of the primary, loggi and the temperature 
of the secondary, 7^2, which are poorly constrained by 
the standards of single white dwarf spectra of compa- 
rable quality. Presumably the simultaneous fit of two 
spectra introduces an element of degeneracy not present 
for single white dwarfs. 

3.4.3. SDSS and UVOT spectral energy distribution 

At optical wavelengths the cool, low-mass primary 
dominates the observed flux of SDSSJ1257+5428. How- 
ever, given the substantial difference in effective temper- 
atures, the flux contribution of the two white dwarfs is 
expected to reverse in the ultraviolet, which motivated 
our Swift TOO observations. These show an increasing 
flux towards shorter wavelengths (Fig.[9])- We have fit- 
ted the energy distribution spanned by the broad-band 
fluxes in the UVOT uvw2, uvm2, uvwl, and b-bands 
and in the SDSS ugriz-bands with composite white 
dwarf models, using the same model grid as in Sect. l3.4TT1 
and 13.4.21 More specifically, we fold all model spectra 
through the spectral response curves of the Swift and 
SDSS bandpasses, and calculate absolute magnitudes 




Fig. 9.— The spectral energy distribution of SDSSJ1257+5428. 
SDSS and UVOT broad-band fluxes are shown as filled circles, 
also plotted is the SDSS spectrum (gray). The best-fit to the 
broad-band fluxes is shown by the solid black line (middle in the 
FUV), with the individual contributions of the cool and hot com- 
ponents l^^dashed_gr^^lme^ ; _A cran^osite model for the param- 
eters of Kulkarni & van Kcrkwiik (2010) is shown in blue (highest 
in the FUV ); our best fit to the phase-dependent WHT spectra 
(Sect. 13.4.21 Fig. in red (lowest in the FUV). 

using the cooling sequences of Bergeron et al.l (H995fJ 
which we extrapolated down to log g = 6.0. In a fi- 
nal step, the co-added set of absolute magnitudes are 
fitted to the Swift and SDSS data. The best fit for- 
mally gives Ti cff = 6620 ± 200 K, loggi = 6.1 ± 0.1, 
T 2 , eff = 12400'± 400 K, and logq 2 = 8.6 ± 0.2 whic h 
are similar to those of Ku lkarni fc van Kerkwijkl (|2010f ). 
but suffer from the same limitation, namely being hard- 
bounded at log<? = 6.0, in our case by the cooling se- 
quences we used. Nevertheless, the UVOT fluxes clearly 
favor a temperature for the secondary at the high end of 
the values returned by our spectroscopic fits. 

Figure[9] illustrates our average best-fit to the phase- 
binned spectra, the fit to the UV/optical broad- band 
fluxes, an d a model computed using the p arameters of 
Table 1 in Kulka rni fc van Kcrkwiik (|2010f ) along with 
the Swift and SDSS broad-band fluxes, and the SDSS 
spectrum. The parameters used in the models are those 
summarised in Table [TJ The models diverge in the 
far ultraviolet and will be easily distinguishable with 
HST spectra. The model based upon the mean of our 
phase-binned data is too cool, while that based upon 
IKulkarni fc van Kerkwijkl (f2010h is a little too hot in the 
FUV and too cool at optical wavelengths. The tem- 
perature of the secondary from the disentangled spec- 
tra (9800 K) is obviously much too cool and we plot no 
corresponding SED in Fig. |H1 

In summary, the optical and ultraviolet characteris- 
tics of SDSSJ1257+5428 are consistent with a double- 
degenerate consisting of a cool, extremely low- mass white 
dwarf plus a rather massive, hotter white dwarf. We esti- 
mate the distance to the bina ry to be d ~ 100 pc, la rger 
than the 48 pc estimated by iBadenes et al.l (|2009f ) be- 
cause of the larger size of the primary. However, the 
data currently available are insufficient to fully constrain 
the five free parameters (Ti^fi, log Si) ?2,cff7 log gi and 
wsinz). 

1 An updated 2006 version is availabl e at 
http: //www. astro. umontreal.ca/'bergeron/CoolingModels 
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Fig. 10. — From top to bottom the symbols with error 
bars show the temperature and gravity of the low-mass pri- 
mary based upon (i) ou r phas e-resolved spectra (triangle) , (ii) 
Kulkarni & van Kerkwiik (2010)'s equivalent result (circle), and 
(iii) our fit to the disentangled spectra (square). The lines show 
evolutionary mode ls for different mass helium white dwarfs from 
IPanei ct al. (2007|). The cooling ages in Gyr are labelled on the 
0.1604 Mq model which is slowed by hydrogen fusion. The arrow 
in the lower right indicates the lower limit on gravity based upon 
assuming that the secondary contributes at least 25% of the flux. 



4. DISCUSSION 
4.1. Present state of SDSSJ1257-, 

The precise nature of the low mass primary star in 
SDSSJ1257+5428 presents us with a puzzle. As Fig. [TO] 
shows, the fits to the disentangl ed spectra are con sistent 
with a mass of around 0.2 M© (jPanei et al.ll2007t ). The 
cooling age in this case is ~ 1 Gyr which just about allows 
the secondar y star (cooling age 1 to 2 Gyr, for CO an d 
ONe models (|Bergeron et aUll995t lAlthaus et aLll2007j )) 
to have a longer cooling age than the primary, as we ex- 
pect, since the massive white dwarf presumably formed 
before the helium white dwarf. In contrast, our phase- 
resolved fits, which fit the data directly , as w ell as the 
similar fits of iKulkarni fc van Kerkwijkl (|2010f ). suggest 
a much lower gravity and a mass < 0.16 M©. In fact, 
as both Figs [5] and [TO] show, the gravity of the pri- 
mary is very poorly constrained by the spectroscopic 
fits, ranging from logg — 5 to 6.6 (and even then the 
lower bound is set by the model grid). Crucially, this 
range allows the primary star to sit on either side of 
the division at ~ 0.18 M© below which residual hydro- 
gen on the surface of helium white dwarfs can maintain 
fusion and greatly slow their cooling whereas above, hy- 
drogen shell flas hes can remove the hyd r ogen and allow 
normal cooling ([Kippenhahn et anil9 68: Webbink 119751 : 
Driebe et aL|ll998tlSarna et alJ | 2000HAlthaus et al lboOll : 
Nelson et al.ll200l IPanei et al.ll2007t ). As Fig. [TO] shows, 
if the mass of the primary is below this limit, then in or- 
der to keep the cooling age of the primary short enough 
to allow the secondary to be as hot as it is, the gravity 
must lie at the lower end of the range with logg ~ 5. 

Such a low gravity is hard to believe because, for rea- 
sonable assumptions about the masses of the two stars, 
the primary becomes so large that its drowns out the 
light from the secondary. Quantitatively, the ratio r of 
the flux from the secondary divided by that from the 



primary, at wavelength A is given by 

R 2 2 F(X,T 2 ) 
r RfF(X, Ti) ' 



(2) 



where F is the astrophysical flux density (flux per unit 
wavelength per unit area of the star). Using R 2 oc M/g 
and applying the observed 3:1 maximum ratio (from the 
disentangling light ratio limits) gives therefore 



M 29l F(X,T 2 ) 1 
M 1 g 2 F(X,T 1 ) > 3" 



(3) 



This then leads to a lower limit upon the gravity of the 
primary: 



1 MxgiF{X,Ti) 
91 3 M 2 F(X,T 2 ) 



(4) 



Assuming M 2 = 1.0 M©, log52 = 8.64, Mi = 0.15 (a 
lower limit to give the least stringent constraint on gi), 
Ti = 6300, T 2 = 14000 K (an upper limit, again to give 
the weakest constraint on g\), and taking the ratio of 
model atmospheres at A = 483 nm (g-band central wave- 
length), we find log<?i > 6.3. We plot this lower limit on 
the right-hand side of Fig. [TO] This limit suggests that 
the surface gravity of the primary star is at the upper 
end of the range returned by the spectroscopic fits. If 
the primary has a mass below 0.18 M©, as the spectro- 
scopic fits seem to favour, then the evolutionary models 
suggest that its cooling age is of order 10 Gyr. 

We are left with a choice between two unpalatable 
alternatives: either the primary has a mass above > 
0.18M©, as suggested by the fits to the disentangled 
spectra but in contradiction with both our direct fits 
and those of IKulkarni &: van Kerkwilkl ([2010), or it has 
a mass below this limit and is therefore extremely old, 
making the youth of the more massive secondary impos- 
sible to understand, since it should have formed first. 

We do not have a satisfactory resolution of this prob- 
lem, but marginally prefer the higher mass solution 
(Mi ~ 0.2 M©) for the primary to avoid the relative age 
problem. This supposes that the evolutionary models 
are correct; further work to fix Mi is encouraged. If 
we are right that the primary mass in SDSSJ1257+5428 
is 0.2 M©, then the mass of the secondary is > 1 M© 
(Fig. [8]). This puts the secondary close to the dividing 
line between CO and ONe white dwarfs which is thought 
to lie at around 1.05 M© for solar metallicity progeni- 
tors, corresp onding to an initial mass of about 5.5 M© 
(|Meng et al.l[200cl . 

We suspect that the discrepancies in the best-fit pa- 
rameters between the various fits result from difficul- 
ties in establishing reliable continua given the extremely 
broad wings from the high mass white dwarf and, in our 
case at least, wavelength dependent slit losses during the 
rather poor conditions of our observing run. Differences 
between the model atmospheres and data are visible in 
both Figs 2] and [7] As a result, our best-fit parameters 
suffer from some systematic errors, which are only com- 
pounded by the degeneracy from fitting two spectra. The 
independent evidence provided by the UV-optical SED 
favors a high temperature for th e secondary as found by 
Kulk arni fe; van Kerkwiik] ([2010), suggesting that their 
data are less affected by continuum problems than ours, 



although their low gravity for the primary possibly in- 
dicates that they are not entirely immune to the same 
problems. Further observations are needed to clarify this 
issue. 

4.2. Anomalous broadening 

Although there has been debate over the first phase of 
mass transfer in the formation of double white dwarfs, 
the second phase, involving mass transfer onto the white 
dwarf that forms first, is almost universally ass umed 
to lead to the formation of a common envelope (jHanl 
Il998t iNelemans et alJfeOOlbt) . However, the absence of 
a narrow core in Ha of the secondary white dwarf in 
SDSSJ1257+5428 is suggestive of rapid rotation with 
wsini > 500 km s -1 , which throws doubt upon this as- 
sumption since we estimate from consideration of the 
moment of inertia of the white dwarf and the angu- 
lar momentum of the accreted material that it should 
have accreted > 0.02 M Q to attain this rotation rate. 
This would require a much more prolonged period of 
accretion than allowed by a common envelope, given a 
maximum accret ion rate onto a wh ite dwarf of a few 
xlO _7 M yr" 1 (jNomoto et alJl2007t ). This is reminis- 
cent of evolutionary scenarios invo lving stable secon d 
phases of mass transfer discussed bv lSarna et al.l fl996). 
A possible near-relative of SDSSJ1257+5428 is the star 
HD 49798 in which a ma ssive, rapidly spinning white 
dwarf orbits an sdO star (jMereg hetti et al.l[200l . Un- 
certainty remains however over the cause of the broad 
profiles. There is no sign of Zeeman spli t ting i n Ha 
(Fig. [5]), but as iKulkarni fc van Kerkwiikl (|2010ft sug- 
gest, a magnetic field could make the core more diffi- 
cult to see, and, if suitably distributed in strength across 
the surface, would not necessarily lead to obvious split- 
ting. In a dual rotation-plus-magnetic-field model, the 
field could also act to increase the lever arm for accre- 
tion and hence reduce the amount of material needed 
to spin the white dwarf up. Alternatively, if the atmo- 
sphere is helium-dominated, we will see deeper into the 
star which may allow Stark broadening to wipe out the 
Ha core. Finally, the temperature of the secondary star 
is such that it could lie within the ZZ Ceti pulsat ional 
instability strip. It is observed ()Koester et al.lll9"98l) that 
ZZ Ceti stars show weak cores in H a, probably the re- 
sult o f motion in their photospheres ([Koester fc Kompal 
[2003). It would be worth observing SDSSJ1257+5428 
photometrically for pulsations. 

4.3. Future evolution 

Close pairs of white dwarfs are common within our 
Galaxy, and their future evolution has been the subject 
of much discussion, not least as potential progenitors of 



Type la supernovae. For SDSSJ1257+5428, the first step 
is clear: gravitational radiation will reduce the orbital pe- 
riod to the point at which mass transfer starts. The or- 
bital period at this point will be about 1.5 minutes which 
will be reached in ~ 3 x 10 9 years. The outcome of mass 
transfer is uncertain and our constraints allow any of 
three possibilities which are (i) explosion as a Type la su- 
pernova, (ii) accretion-induced collapse to a millisecond 
pulsar and (iii) survival of the onset of mass transfer as 
a hydrogen-deficient ultra-compact binary. We slightly 
favor the third possibility on the grounds that the mass 
ratio q = Mi / M 2 re 0.2 is extreme enough to allow stabl e 
mass transfer (jNelemans et al.l2001aHMarsh et al.l2004ft . 
while we suspect the total mass to be less than the Chan- 
drasekhar mass. SDSSJ1257+5428 is therefore likely to 
become a semi-detached, accreting double white dwarf 
(AM CVn star). The secondary star will become the ac- 
cretor. Since it is massive for a white dwarf, and as long 
as it is not an ONe white dwarf, the system is a good can- 
didate prog enitor of the sub-lum inous ".la" supernovae 
dis cussed bvlBildsten et all (|2007l ) and possibly observed 
bv lKasliwal et al.l (j2010fl 

5. CONCLUSION 

We find that the putative white dwarf-black- 
hole/neutron star binary, SDSSJ1257+5428, is a double 
white dwarf. SDSSJ1257+5428 is composed of a very 
low mass ~ 0.2 M Q white dwarf together with an ex- 
tremely massive (> 1 M Q ) white dwarf. As long as the 
massive white dwarf avoids accretion-induced collapse or 
explosion, SDSSJ1257+5428 will evolve into a hydrogen- 
deficient accreting binary star, but may later explode as 
a sub-luminous Type la. The massive white dwarf shows 
signs of rapid rotation (vsini > 500 kms -1 ) which sug- 
gests that the most recent phase of mass transfer might 
not have involved a common envelope, contrary to cur- 
rent models of double white dwarf populations. Some 
inconsistencies in the parameters of the two white dwarfs 
remain that are probably caused by difficulties in fitting 
their blended spectra; further observations are required 
to clarify these. 
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tained by the Centre Donnees astronomiques de Stras- 
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the authors. 
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